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DOUBLE RESONANCES AND SPECTRAL SCALING IN THE WEAK TURBULENCE 
THEORY OF ROTATING AND STRATIFIED TURBULENCE 

ROBERT RUBINSTEIN* 

Abstract. In rotating turbulence, stably stratified turbulence, and in rotating stratified turbulence, 
heuristic arguments concerning the turbulent time scale suggest that the inertial range energy spectrum 
scales as k ~ 2 . Prom the viewpoint of weak turbulence theory, there are three possibilities which might 
invalidate these arguments: four- wave interactions could dominate three- wave interactions leading to a 
modified inertial range energy balance, double resonances could alter the time scale, and the energy flux 
integral might not converge. It is shown that although double resonances exist in all of these problems, 
they do not influence overall energy transfer. However, the resonance conditions cause the flux integral for 
rotating turbulence to diverge logarithmically when evaluated for a k~ 2 energy spectrum; therefore, this 
spectrum requires logarithmic corrections. Finally, the role of four-wave interactions is briefly discussed. 

Key words, rotating turbulence, stratified turbulence, weak turbulence 

Subject classification. Physical Sciences 

1. Introduction. Weak turbulence theory [1, 2, 3] provides a closure for problems of interacting disper- 
sive waves. As developed by Zakharov and collaborators [3], it provides a means to calculate Kolmogorov-like 
spectra carrying constant fluxes of invariant quantities, to decide whether these spectra are independent of 
the large-scale excitation and the dissipative mechanism (the problem of locality ), to decide whether energy, 
or other invariants, are carried from large to small scales or from small to large scales (direct or inverse 
cascades), and to assess the stability of these spectra to various perturbations. A two-point closure closely 
related to weak turbulence theory has been applied by Cambon and co- investigators to the numerical simu- 
lation of rotating [4] and of stratified [5] turbulence. 

An important ingredient in this theory is the resonance operator 

/ oo 

dr e i^(k,P,q) 5(k _ p _ q) 

-oo 

where h is defined in terms of the dispersion relation uj = u>(k) of the waves by 
(1.2) h(k , p, q) = cj(k) ± w(p) ± w(q) 

Eq. (1.1) has been written for a quadratically nonlinear theory, in which the lowest order perturbation theory 
[2, 3] leads to three- wave interactions. Note that each choice of signs in Eq. (1.2) leads to a different operator. 

The resonance operator defined by Eq. (1.1) depends on the geometric properties of the resonance surface 
h = 0. The most important question about the resonance surface is whether it contains any (real) points 
at all. If /i(k, p, q) never vanishes when k = p + q, then the integrand in Eq. (1.1) oscillates strongly and 
the resonance operator vanishes. If all three-wave resonance operators vanish, it is necessary to proceed to 
the next order in perturbation theory, which generates a theory with a nonlinearity of higher order than 
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the orginal equations of motion; for quadratically nonlinear problems, the result is a theory with a cubic 
nonlinearity. In such a theory of four-wave interactions, the resonance surface is never empty [3]. 

Assuming that the resonance surface is nonempty, the next question is whether it is nonsingular. If it 
is, then elementary partition of unity arguments [6] show that 

/ oo 

dr e iT M k .p>q) = S (h) 

-OO 

Consequently, if the dispersion relation is homogeneous of degree a, 

(1-4) w(Ak) = A Q w(k) 

then 

(1-5) 7£(Ak, Ap, Aq) = X~ d a TZ{k, p, q) 

where d is the dimension of space. Eq. (1.4) states that the turbulent time scale is of order k~ a . For problems 
with three-wave interactions, including problems with hydrodynamic nonlinearity, a simple scaling argument 
given by Orszag [7] suggests that when the energy flux is constant, the spectrum scales as 

(1-6) E(k)~k~ 2+a/2 

A complication which these arguments ignore is the possibility that all three-wave resonance operators 
vanish. As noted earlier, turbulent energy transfer is then dominated by four-wave interactions and the 
inertial range energy balance which leads to Eq. (1.6) must be modified. For later reference, we note the 
energy spectrum corresponding to four- wave interactions 

(1.7) E(k) ~ k ~ 7 ' 3 +“/3 

Another complication which is ignored by the heuristic arguments is the existence of double resonances , 
which arise whenever the resonance surface is singular. The arguments leading to Eq. (1.3) are no longer 
valid, and the theory leads to a more complex resonance operator and possibly to a different time scale. 
The classic case of double resonance is nonlinear sound waves [8, 9]. In two space dimensions, the double 
resonances change the time scale of the problem completely, and the scaling suggested by Eq. (1.6) does not 
apply [3]. This problem will be discussed in more detail in the next section. 

The dimensional arguments leading to Eq. (1.6) also assume that no time scale except the linear dispersive 
time scale is relevant, and that no length scale except the wavelength ~ A: -1 is relevant; constancy of the 
energy flux then suffices to determine the spectral scaling. In particular, neither the production nor the 
dissipation mechanism influences the spectral scaling. This locality condition noted earlier can be expressed 
analytically as the convergence of an integral representing the energy flux when this flux integral is evaluated 
for an infinite inertial range given by Eq. (1.6). Divergence at large or small scales indicates the relevance of 
the production or dissipation mechanism respectively, and alters the spectral scaling. Locality depends on 
the geometry of the resonance surface, since the flux integral might diverge along singular curves or points 
of the surface. 

This paper considers three related problems of weak turbulence theory: rotating turbulence, stratified 
turbulence, and the coupled problem of rotating, stratified turbulence. In all three problems, the dispersion 
relation is homogeneous of degree zero: 

(1-8) w(Ak) = w(k) 
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In the notation of Eq. (1.4), a = 0; consequently Eq. (1.6) predicts 

(1.9) E(k) ~ k~ 2 

The k~ 2 spectrum has been already been proposed for rotating turbulence by Zhou [10] and for stratified 
turbulence by Herring [11] using arguments similar to those of Orszag [7] noted earlier. Although all three 
problems contain a preferred direction, and must exhibit spectral anisotropy, we consider in Eq. (1.9) the 
energy averaged over shells in wavevector space. 

The question which we pose is whether the scaling Eq. (1.9) actually is the prediction of weak turbulence 
theory for these problems. In order that this be so, it is necessary to 

1. evaluate the effect of any singularities of the resonance surface 

2. demonstrate the convergence of the flux integral. 

3. demonstrate the applicability of three-wave theory 

It will be shown that although double resonances exist in all of these problems, they do not influence 
energy transfer; consequently, the double resonances do not alter the energy spectrum. But evaluation 
of the energy flux integral for rotating turbulence shows that the k~ 2 spectrum is only marginally local, 
and that logarithmic corrections analogous to those derived by Kraichnan [12] for the enstrophy range in 
two-dimensional turbulence are necessary. Accordingly, instead of Eq. (1.9), 

(1.10) E(k) ~ fc- 2 [log(fcLo)]- 1/2 

where Lq is the integral scale. Finally, the role of four- wave interactions in these problems is briefly addressed. 
These problems admit three- wave interactions, except for one special case of coupled rotating, buoyant 
turbulence. 

Since the k~ 2 spectrum also occurs in Burgers turbulence, a connection is sometimes made between the 
occurrence of this spectrum and the presence of wave breaking [13]. The physical origin of the k~ 2 spectrum 
in weak turbulence theory is certainly unrelated to wave breaking, which is an effect of strong nonlinearity. 
In this connection, it is interesting that Kraichnan’s derivation [14] of the spectrum of Burgers turbulence 
does not invoke wave- breaking, although it is also a strong turbulence result. 

The analysis of this paper will be limited to wave interactions alone: it will only consider the turbulence 
of waves satisfying the dispersion relation appropriate to each problem. This is certainly a limitation in the 
cases of stratified turbulence and the coupled problem, in which interactions with non-oscillatory geostrophic 
modes [13] dominate the dynamics [15]. Since there is no three- wave mechanism which generates these modes 
from wave interactions [13], this analysis appears to be self-consistent. However, the question is very much 
open whether four- wave interactions cannot generate geostrophic modes [16], an interaction which would 
have considerable importance in geophysical applications. 

2. Double resonances in weak turbulence theory. The classic case of double resonances in weak 
turbulence theory occurs for nonlinear sound waves [8, 9], where 

(2.1) w(k) = c|k| 

Double resonances occur if 

Px _|_ q 

P Q 

?y±?y= o 

P Q 

(2.2) — ± — - 0 

V Q 
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These conditions are equivalent to 


(2.3) 


±E = El = Ey _ p± 

Q Qx Qy Qz 


and can be satisfied for any k if p, q and k are collinear. 

In three dimensions, the double resonance induces a logarithmic time scale which requires a re-ordering 
of the perturbation theory [8] on which weak turbulence theory is based, but it does not to alter the scaling 
relation Eq. (1.5) or the spectral scaling, E(k) ~ A; -3 / 2 which Eq. (1.6) suggests since u ~ Jfe 1 . 

In two dimensions, however, the situation is quite different: the resonance operator proves to scale with 
time as t 1 / 2 , and the theory must be reformulated completely [9]. Thus, if the t 1 ? 2 time dependence is 
absorbed into the definition of a new slow time scale, the resonance operator would scale as 1Z ~ (k/c) 1 ^ 2 
instead of as 1Z ~ k/c, leading to E(k) k* 7 / 4 . 

The problem of nonlinear sound waves is exceptional in several respects [3], and the correct spectral 
scaling is now considered to depend on higher order terms in the dispersion relation [17]. Nevertheless, 
this example illustrates how double resonances can modify the turbulent time scales and consequently alter 
the predictions of weak turbulence theory. We next determine whether singularities exist in the resonance 
surface for the three problems of this paper. 


3. Resonance surface singularities. The three cases of rotating, stratified, and rotating stratified 
turbulence will be considered in turn. Although rotation and stratification are formally special cases of the 
coupled problem, they prove to have special features requiring separate analysis. 


3.1. Rotating turbulence. The dispersion relation of inertial waves in rotating turbulence is 

(3.1) w(k) = ±2Q~ 

k 

where Q is the rotation rate. The resonance surface is therefore 


(3.2) h(\t, p, q) = ^± — ± — = 0 

k p q 

Note that the resonance surface is nonempty and three-wave interactions are possible in this problem. 

We must consider the possible singularities of the surface h = 0 when k, p, q satisfy the triangle condition 
k = p + q. Singularities, or double resonances, occur if dh/dp = 0, or 


(3.3) 

when h = 0. Eq. (3.3) states that 


d ,p z q Zx 

— (Z*±**) = 0 

dp p q 


(3.4) 


PxPz ,QxQz 
p 3 * q 3 

PyPz , QyQz 

p 3 ~ *”^ 3 " 

p * + p2y = + + q y 

p 3 ± q3 


It is understood that the sign in each equality in Eq. (3.4) is the same, hence only two cases are represented. 
But the last equality can only be satisfied if both sides are positive; therefore, the negative sign case in 
Eq. (3.4) does not give any double resonance. 
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Assume first that no component of p vanishes. Rearranging the terms, 


(3.5) 


£ _ 9x9; _ 9y9* _ 9x +9y 

p 3 PxPz PyPz p\ + Py 


therefore 

(3.6) 


9x _ 9y _ 9^Px 
Px Py P 3 Qz 


and 

(3.7) 

consequently 


£ pi 
p 6 qI 


(3.8) 


Qx Q z 

Px Py Pz 


The homogeneity of the dispersion relation means that the first set of equalities can only be satisfied if 


Note the role of the homogeneity of this problem uj ~ k° in forcing this conclusion: the result Eq. (3.9) does 
not apply to sound waves, because u> ~ k l . 

We conclude that the conditions for double resonance can only be satisfied if either q — p or q = - p. In 
the first case in which k = 2p. the wave vectors are collinear, and the dispersion relation requires ±1±1±1 = 0, 
which is impossible. The second case in which k = 0 pertains to a degenerate mode which we do not consider. 
It follows that double resonances can only occur if some component of p vanishes. 

Next, consider these special solutions of Eq. (3.4). There are two obvious solutions which satisfy the 
resonance conditions. First, if p z — q z = 0 and p = q, the resonance conditions are trivially satisfied, since 
also k z = 0. Under these conditions, p = q implies that p lies on the line 


(3.10) 2p ■ k = k 2 ,p z = 0 

If k z = 0, the resonance conditions require that either p z = q z = 0, or that 2p • k = k 2 . These planes 
intersect in the line defined by Eq. (3.10), which is therefore a double line of the resonance surface. 

A second obvious solution of Eq. (3.4) is p x = = 0, where (p^) 2 — p 2 -hp 2 . In this case, the resonance 

condition cannot be satisfied if k is a third vertical vector, because, as noted above, the resonance condition 
and collinearity would imply ±1±1±1 =0. The only possibility is that k z = 0. But then the triad condition 
k = p + q forces k = 0. Therefore, we again are led to resonance with a degenerate mode, a case which can 
be ignored. 

The remaining instances of vanishing components do not lead to double resonances not already noted: 
the elementary case-by-case verification appears in Appendix I. Accordingly, we conclude that the resonance 
surface has exactly one singularity, the double line defined by Eq. (3.10). 

3.2. Waves in stratified turbulence. Internal waves in stratified turbulence satisfy the dispersion 
relation 

k 1 

(3.11) w(k ) = ±N — 

fa 
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Three-wave interactions are possible for this problem. 
For a double resonance, 


(3.12) 

and h — w(p) ± u;(q) ± u;(k) = 0. As in Eq. 
first two equalities in Eq. (3.12), 


PxP'j = ± Mi 

p^p^- Q 3 Q^ 

PvP\ = , QyQz 
p3p± q 3 q ± 

PzP X _ q z q x 

pi ± q 3 

(3.4), only two cases are represented by Eq. (3.12). Prom the 


(3.13) 

Using the third relation of Eq. (3.12), 


Px = Py = ± P^P^<£ 

9* q y q 3 q L p 2 


± ^_ptsl = +pL'1l'p^Sl - 2<k. 

q 3 q ± P 2 z Pz q 3 q 1 ' Pz V J Pz 


where only the positive sign applies in the last equality. Since 


(3.15) 

the last equality implies also 

(3.16) 

and finally, as in Eq. (3.8), 


Px = Py _ ± p^_ 

qx q y q 1 - 


= = P| (h 

q 3 q ± pi { q ± Pz qlp z 


(3.17) Px _ Py = Pz 

Qx Qy Qz 

Since the homogeneity implies that necessarily q/p — 1, we are again led to the two cases considered in 
the previous section. The dispersion relation again excludes the case of collinearity p = q and k = 2p. If 
q = -p, we obtain a degenerate mode k = 0 and this case is again excluded. 

Next, consider the solutions of Eq. (3.12) when at least one component of p vanishes. First, if p z = q z = 0, 
then also k z = 0, and the resonance condition cannot be satisfied since it requires ±1±1±1 — 0. 

Second, if = q 1 - = 0, the third equation is satisfied but the first two become indeterminate. But it 
is clear that three vertical vectors satisfying the triad condition are always resonant. Such vectors define a 
degenerate line component of the resonance surface; any such component must be singular. 

The remaining cases of vanishing components do not lead to anything new, and are analyzed in Appendix 
II. 


3 . 3 . 


(3.18) 


Waves in rotating stratified turbulence. The dispersion relation is 

W (k) 2 = 4fi 2 (h)2 + JV 2 (^)2 
rC K 

= ( 4fi 2 -7V 2 )(h ) 2 + j V 2 
K 
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Three-wave interactions are possible provided 2 D ■/ N. When 211 = N, two waves can only interact with a 
geostrophic mode [13] and interactions among waves alone must be described by a four-wave theory. 

Double resonances occur provided 


(3.19) 


J_PiP| = ± 1 9*9z 
u p p A uj q q A 

1 PyPz _ ± 1 9t/9z 

P 4 Wg g 4 

1 P ^ 1 ) 2 J. 1 ^( g±)2 

Wp p 4 Uj q q 4 


As in the two previous cases, 

(3.20) 

Using the third relation from Eq. (3.19), 

(3.21) 


Px Py j Qz 

q x q y u q q 4 p 2 z 

1 - 

r, r,4 


consequently, 

(3.22) 


Px _ Py _ 9z ,P\ 2 _ 9zPx 
9i 9s/ p 2 9 1 Pz 9* 


so that the conclusion reached in the two special cases, 


(3.23) 


Px _ Py_ _ Pz _ 

Qz 


again holds. Since the homogeneity of the dispersion relation requires that the factor of proportionality in 
Eq. (3.23) have magnitude one, we again come to the two cases before, neither of which leads to a double 
resonance. 

The solutions of Eq. (3.19) when at least one component of p vanishes include first p z = q z = 0. Then 
also k z = 0 and the resonance conditions cannot be satisfied since they reduce to±l±l±l = 0. The same is 
true of the second obvious solution = q "*■ = 0. Thus, neither of the cases which lead to double resonances 
in the uncoupled problems remains in the coupled problem. 

However, the solution p x = p y = 0, q z = 0 does lead to a double resonance. In this case, the resonance 
condition reduces to 


(3.24) 


4f } 2 p 2 z + N 2 (q ± ) 2 = k 2 { 2n ± AT) 2 


and the choice of the negative sign gives the equation 

(3.25) 0 = k 2 2 (N 2 - 4DiV) + (k 2 x + ^)(4D 2 - 4 SIN) 

since k = ( q x ,q y ,Pz )■ Thus, for every k satisfying Eq. (3.25), the decomposition p = (0, 0, p z ), q = (k x . k y , 0) 
defines a double resonance. Eq. (3.25) is the equation of a cone provided N > ft or N < 4fl; if fl < N < 4D, 
then there are no double resonances at all. The second alternative applies to the case of four-wave resonance 
N = 2fi noted earlier. If we consider the general case in which q z = 0 and look near the point Px = Py — 0 
when k satisfies Eq. (3.25), it is found that the singularity is that of a line component, the same singularity 
that occurs in stratified turbulence. 

The remaining case in which Eq. (3.19) can be satisfied is p x = q x = 0. As in the previous cases, this 
reduces to the case of collinearity, which has already been considered. 
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4. Singularities and the resonance operator. To determine the effect of the multiple resonances, 
it is necessary to consider the form of the resonance operator near the singularity. For example, in the case 
of nonlinear sound waves, all resonances are double because the resonance condition forces colhnearity of the 
three wavevectors. Locally, the resonance operator behaves like the generalized function 

( 4 - ! ) [X?,/] = f dV f dre l( - x '+y^ r f(x,y,z) 

j J — oo 

in which the resonance surface has degenerated into the line x = y = 0. In Eq. (4.1), fix, y. z) is an arbitrary 
test function. Changing to polar coordinates, 

/ OO rOC yOO y2 7T 

dr dz rdr / dOe"^ f{r,6,z) 

-oo */ — oo J 0 JO 

/ oo 

dzf{ 0,0, z) 

-oo 

Consequently, [8] the resonance operator reduces to integration along the line representing the degenerate 
resonance surface. The three distinct integrals in Ref. 8 correspond to the different arrangements of three 
collinear vectors: k > p > q, p > k > q, or p > q > k. 

This same approach will be applied to analyze the effect of double resonances in rotating and stratified 
turbulence: the resonance operator will be represented locally by a phase function with the same singularity 
as the resonance surface. Since the double locus is different in rotating and stratified turbulence, each case 
must be analyzed separately. 

4.1. Rotating turbulence. Unlike the case of nonlinear sound waves, the resonance surface is non- 
degenerate for rotating turbulence, and it is nonsingular away from the double curve. Consequently, one 
contribution to the resonance operator is just £(/i). The problem remains to determine the effect of the 
double curve on the resonance operator. 

Near a double curve, the resonance operator behaves like the generalized function 

( 4 - 3 ) ^ V ' f ^ = J dV f 0 rfre<IWT /(x,3/,2 ) 

in which x = y = 0 is a double line of the surface xy = 0. In polar coordinates, 

roc p OO poo p'2tT 

(4.4) [£>,/]= / dr dz I rdr / d0e ir7 % 8in ( 20 ) T /( r , 0, z) 

Jo J —oo Jo Jo 

To determine the effect of the double curve, it will suffice to ignore the variation along the double curve and 
evaluate 



in the special case / = cxp(- ^ar 2 ). If the limit 


( 4 -6) lim [P*, e - ar2 l 

a—* oo J 

is zero, the double curve makes no contribution to the resonance operator; if it is not zero, then there is an 
additional contribution corresponding to integration over the double curve. 

Integrating over r, 


(4.7) 


[D\e 


p2n pa 

= / do 

Jo Jo 


dr- 


a — irsin(20) 
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The imaginary part of this integral is logarithmic, but it makes no contribution to the energy balance [2]. 
The real part of the integral in Eq. (4.7) can be evaluated as 


(4.8) 


Re[D* , e 


= lim 8 
r— oo 


/ 


dO 


7t/4 

t /4 sin(20) 


arctan(— sin(20)) 
a 


The integral can be evaluated approximately by setting 


(4.9) 


arctan(— sin(20)) 
a 


T sin(2 9) /a 
tt/2 


if | sin(20) |< 7ra/2T 
if | sin(20) |> na/2T. 


Substituting Eq. (4.9) in Eq. (4.8), 

(4.10) [D m ,e~ ar2 ] « 7T + j lim c 27rlog(^) 

This approximation defines an upper bound, since the left side of Eq. (4.9) is smaller than the right side. 

We conclude that the double curve contribution to the resonance operator is logarithmically divergent 
in time. Therefore, for rotating turbulence, the real part of the resonance operator is 

(4.11) Re'll = ^ ^2 <*M k ) ± w(p) ± w(q))<5(k - p - q) + T*6{C) 


where C denotes the double curve Eq. (3.10). The notation T* is used provisionally to represent the 
logarithmically divergent time scale. 

In terms of the multiple-scale perturbation expansion used in weak turbulence theory, the result in 
Eq. (4.11) should be understood to mean that there are two slow time scales in rotating turbulence, 

(4.12) Ti = e 2 t, T 2 = e 2 J log (s)ds 

where e is the expansion parameter [2]. Unlike the otherwise analogous case of sound waves, this logarithmic 
time scale appears in the real part of the resonance operator and therefore could influence the inertial range 
energy balance. 

However, there is an important difference between the present case and sound waves: whereas in the 
problem of nonlinear sound waves, double resonances occur for all possible k, in rotation, they occur only for 
horizontal vectors. The weak logarithmic singularity induced by the double resonance will disappear upon 
integration over all vectors k in the flux integral. Moreover, the double resonance occurs for interactions 
which do not transfer any energy into the mode k [18]. We conclude that the double curve makes no 
contribution to the energy flux. 


4.2. Waves in stratified turbulence and rotating stratified turbulence. In both problems, the 
Taylor series of the resonance surface near the singularity begins with the quadratic term p \ + p^ which is 
satisfied by the line p x — Py — 0. This is also the singularity found in the case of nonlinear sound waves. The 
local behavior of the resonance operator is therefore given by Eq. (4.1) and the conclusions of Newell and 
Aucoin [8] apply: the double resonance causes a logarithmic time scale, but does not modify the spectral 
scaling. 

But, as in the case of rotation, the logarithmic time-scale arises only for a subset of vectors k satisfying 
at least one equation: vertical vectors in the case of stratified turbulence, and the condition Eq. (3.25) in 
the coupled problem. In forming the flux integral, the logarithmic contributions to the time-scale will be 
suppressed by integration over all vectors k. 

We conclude that the double resonances are are irrelevant to energy transfer and spectral scaling in all 
three problems. 
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5. Convergence of the flux integral. In rotating turbulence, it is convenient to express the flux 
integral in terms of the Craya-Herring basis 


(5.1) 


eW(k) = kxfi/ | kxfi | 
e (2) (k) = kx(kxfi)/ | kx(kxfi) 


or the equivalent basis of Cambon and Jacquin [4] and the corresponding tensors 


%=e^e?-ey>e\ 

th = + 

C 2 _ J1)J1) 


(!)J2) 


e ) e -' 


t *. — e ^v i; - p ( 2 ) p (2) 

Sjj e i c j e , e j 

In Eq. (5.1), 17 = 2(0,0,17) is twice the angular velocity. Note that 


(5.2) 


(5.3) 




Pij( k) = 6ij - kikjk 


-2 


The linear response function for rotating turbulence is the solution of 
(5-4) Gij + PipSlpqGgj = S(t - s)Sij 

which can be written 


Gij(k,t,s) = {cos(2fifc 2 (f — s))Pij(k) 

( 5 - 5 ) + sin(2fifc 2 (f - s)/A:)^° (k)}H(i - s) 

where H is the unit step function. At this level of approximation, the fluctuation- dissipation theorem, 

( 5 - 6 ) r) = Gj m (k, r)Q mi (k) + G jm (k, -r)Q mi (k) 

applies, where r = t — s is time difference. 

It must be expected [4] that the single-time correlation function is a sum over all of the f tensors of 
Eq. (5.2) 

(5- 7 ) Qii(k) = ^C? p (k)^(k) 

But substituting Eq. (5.7) in Eq. (5.6), we find that the symmetry 

Qij(k.,r) = Qji(k, — r) 

leads instead to the simpler expression 

( 5 * 8 ) QijiKr) = sin(2nfc,r)Q(k)^.(k) + cos(2n**r)Q(k)£?.(k) 

in which the single-time correlation function is tensorially isotropic. For the purpose of computing the flux 
integral, we will assume isotropy in k as well. Although this description of the correlation function is certainly 
not realistic, the only use we make of it is substitution in the energy flux balance; an isotropic form is adequate 
for this purpose. In order to evaluate the correlation function, it is necessary to solve closure equations; a 
complete discussion is given by Cambon and Jacquin [4], although the analytical solution remains unknown. 
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The energy flux integral, representing the net energy flux across the wavenumber k' is 
£ = { [ dk f dpdq - f dk f <fpdq}<S(k - p - q) f dr x 

J k<k' ''P,q>k' J k>k' Jo 

{Pim'n (k)-P m'rs (p)Qir (k, ^") Qns (Qj (p? ^*) 

+ 7 ? imn(k)-P mr s(q)Qtr(k, T)Q n5 (p, T)G mTn ' (q, t) 

(5.9) — PimnityPi'rsO^Qmr (Pi ' r )Qns(q> T)Ga f (k, t)} 

for three-wave interactions, where 

(5.10) Pimn (k) = k m Pi n {k) + k n Pim (k) 

Locality of the inertial range means that this integral is finite when the single-time energy spectrum scaling 
as Q(k) ~ fc -4 is substituted. In the absence of resonances, the convergence of the flux integral in this case 
is well-known: [7] it is a special case of the convergence of the integral in Eq. (5.9) when Q(k) ~ and 

1 < (5 < 3. However, the resonance conditions require that this condition be checked independently. 

5.1. Rotating turbulence. The proposed spectrum falls off rapidly enough at small scales that the 
only possible divergence occurs when q 0 and p -> k. The first question is whether the cancellation of 
singularities in the absence of resonance conditions survives if the more complex tensor forms of Eq. (5.8) 
are substituted. The verification is elementary: in this limit, the singularities in Eq. (5.9) are multiplied by 

(5.11) a«»(k)p mr .(k)eS?(k)e«)(o) = knk.PM^^o) 


and 

(5.i2) p, m „(k)p irs (k)^ r (k)^(o) = k n k s p im ( k)eE>e£?(o) 


In Eqs. (5.11) and (5.12), i and j can each equal 0 or 3. When j = 0, both terms vanish, and when j = 3, 
both are proportional to k 2 . Consequently, the terms of order q° cancel. Then the integrand in Eq. (5.9) 
is of order q. But terms of this order vanish on integration by symmetry. Accordingly, as in the case of 
isotropic turbulence, the integrals behave like 


(5.13) 



near q = 0 in the absence of resonance conditions. This result depends on the odd parity of the integrand 
under the substitution q — ► — q; this parity is preserved by the resonance conditions which are homogeneous 
of degree zero. 

The resonance conditions imply that q does not approach zero in an arbitrary fashion. Instead, the 
following cases occur as q — > 0 and p — ► k: 


(5.14) 


cj(k) +u>(p) = 2k,u;(q) = ±2cj(k) 
u;(k) - u>( p) = O.cj(q) = 0 


In the first case, q varies on a cone q z = Xq , and in the second case, q is on a surface tangent to the plane 
q z = 0. In both cases, the element of area on the surface can be taken proportional to qdq. The asymptotics 
indicated in Eq. (5.13) are replaced by 


(5.15) 
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and the flux integral is logarithmically divergent when (3 = 2. 

It is also necessary to verify that there is no additional divergence contributed by the double curve 
Eq. (3.10), but this conclusion is immediate because the double curve does not meet the singularity at q = 0. 

Thus, the flux integral is logarithmically divergent for the k~ 2 spectrum. It is known [12] that in this case, 
logarithmic corrections to the spectrum are necessary: instead of Eq. (1.9), the corrected scaling Eq. (1.10) 
applies. 

5.2. Waves in stratified turbulence. It was suggested earlier that the theory of waves in stratified 
turbulence decoupled from geostrophic motion may not occur in the actual physical system if four-wave 
interactions can generate geostrophic modes from wave modes. Accordingly, a model problem will be inves- 
tigated in which the nonlinearity of rotating turbulence is retained, but the dispersion relation is replaced 
by the dispersion relation of stratified turbulence. 

The conditions in Eq. (5.14) apply to this model problem. In the first case, q approaches zero along a 
cone q x = A g, but in the second case, q must approach zero tangent to the line q x = q y = 0. Approach 
to q = 0 along a line would be extremely singular, but the tangency forces the ratio q ± /q to be nonzero, 
contrary to the assumption that q^/q = 0. Thus, the flux integral is only logarithmically divergent, as in 
the case of rotation. The log-corrected spectrum Eq. (1.10) will apply. 

Analogous arguments apply to the coupled problem. However, we stress that a complete treatment of 
these problems requires an account of the effect of coupling between waves and geostrophic modes. 

6. The role of higher order resonances. It was noted earlier that the energy spectrum 
(6-1) E(k) = Cyffek- 2 

where the frequency f is 2 Cl for rotation and TV for stratification has already been proposed for rotating 
turbulence, [10] and for stratified turbulence [11]. In the coupled problem, Eq. (5.5) has also been proposed 
with a suitable function / = f{Sl,N) [19]. It was noted earlier that in all cases, the heuristic argument is 
similar to Orszag’s [7]: if the time scale is homogeneous of degree zero, simple closure arguments assuming 
a constant energy flux imply that the energy spectrum scales as k~ 2 . 

Weak turbulence theory suggests that such arguments are incomplete. In the isotropic problem of 
Langmuir turbulence, u ~ A; 0 , but the dispersion relation does not permit resonant three- wave interactions 
[3]. In the absence of three- wave resonances, four- wave processes must be considered. This change modifies 
the inertial range energy balance. The flux integral Eq. (5.9) for three- wave interactions can be written for 
scaling purposes as 

( 6 -2) e ~ J(dk) 2 P(k)eP(k)Q(k) 2 

where P ~ A: represents the nonlinear coupling, 0 is the time scale, and the number of wavevector integrals 
is represented symbolically. For four-wave interactions, instead 

(6-3) e~ J (dk) 3 P(k)(QP(k)) 3 Q(k) 3 

If 0 ~ A: 0 , then Eq. (6.2) predicts E(k) ~ k~ 2 but Eq. (6.3) predicts E(k) ~ k~ 7 / 3 . More generally, the 
spectral scalings of Eqs. (1.6) and (1.7) follow from the flux integrals Eqs. (6.2) and (6.3). 

Yakhot [20] has proposed that four-wave resonances play a crucial role in rotating turbulence, especially 
in the small aspect-ratio systems investigated by Smith et al. [21]. It is possible to assess the role of four- 
wave interactions heuristically from the scaling forms of the energy balance Eqs. (6.2) and (6.3). Suppose 
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that the modal amplitudes are determined by three-wave interactions as proposed above. These amplitudes 
can be substituted in the four-wave expression, although different sets of modes resonate in each case. The 
resulting four-wave correction to the energy flux is of order fc 2 e/fl 3 . Since e 1 / 3 ^ 2 / 3 /!! << 1 is the condition 
for weak turbulence theory to apply, the corrections to energy flux due to four-wave interactions would seem 
to be small. 

But this discussion ignores the constants in front of the two flux integrals. In cases in which three-wave 
interactions are depleted, say by a combination of small aspect-ratio and a very large horizontal viscosity, 
the four-wave contribution could dominate the three-wave contribution. 

In stratified turbulence, the resonance condition again permits three-wave resonances, and it has been 
shown that the double resonances are irrelevant to energy transfer. But in this problem, the importance 
of interactions with geostrophic modes [5] and the comparatively secondary role of pure wave turbulence, 
suggests that a more detailed analysis is required. The turbulence of nearly horizontal wave modes has 
recently been analyzed by Caillol and Zeitlin [22] on the basis of weak turbulence theory. 

The same considerations apply to the coupled problem. But in this case, the dispersion relation Eq. (3.18) 
degenerates when N = 2fl to the isotropic form u = N which excludes three- wave resonances. Wave 
turbulence in this case must be analyzed in terms of four- wave interactions, but since any four wavevectors 
can form a resonant quartet, it is not possible to derive a spectral scaling from the flux integral in the weak 
turbulence approximation without making further assumptions. In view of the greater dynamic importance 
of the interactions of waves with geostrophic modes, the role of four-wave interactions in rotating stratified 
turbulence remains uncertain. 

Appendix I. If p x — 0, then either q x = 0 or q z = 0. If p x — Qx = 0, then if p y ^ 0, dividing the 
last condition in Eq. (3.4) by the second, p z /p y = q z /q y . This reduces to the case of proportionality of the 
wavevectors, which we have seen does not yield a double resonance. If instead, p y = 0, then the third of 
Eq. (3.4) forces also q y = 0, and we have two vertical vectors. This case has also been shown not to lead to 
double resonances. 

Suppose next that p x = 0 and q z = 0. Then either p y - 0 or p z = 0. In the first case, p x — Py 0 
implies q x = q y = 0, which does not lead to a double resonance. In the second case, p z = q z = 0, and this 
case has already been considered. 

Other cases correspond to interchanging x and y. The only double resonance is therefore the double line 
Eq. (3.10). 

Appendix II. If p x = 0 and q x = 0, then p y = ±p ± and q y = ±q ± . Multiply the second relation of 
Eq. (3.12) by p^/p 1 - on the left side and by q ± /q- L on the right side, then use the third relation to conclude 
that Pz/Py = Qz/qy This is the case of collinearity, which we know cannot satisfy the resonance condition. 

If Px = 0 and q z = 0 then either p z = 0 or p y = 0. In the first case, p z = q z = 0 which has already been 
analyzed. In the second case, we have p = (0,0,p 2 ),q = iq x , Qy 0) and then p,q, and k cannot resonate. 

The remaining cases again correspond to interchanging x and y. The point singularity found in the main 
text is therefore the only multiple resonance. 
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